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1. 


Introduction, 


In  practice,  long  life  items  are  subjected  to  larger  than  normal 
stresses  in  order  to  obtain  failure  data  in  a short  amount  of  test  time. 

Such  tests  are  called  accelerated  or  overstress  life  tests,  and  the  goal 
is  to  mahe  inferences  about  the  life  distribution  of  the  items  at  the 
normal  stress  levels  using  failure  data  from  accelerated  tests.' 

The  current  approach  to  this  problem  involves  making  assumptions 
about  the  distribution  of  failure  times,  and  about  the  functional  rela- 
tionship between  the  parameters  of  the  failure  distribution  and  the 
3>pplisd  stress.  Such  a relationship  is  known  as  an  acceleration  function 
or  a time  transformation  function;  examples  of  these  are,  the  Power  Law, 
the  Arrhenius  Law,  the  Eyring  Law,  etc. . Another  assumption  commonly  made 
states  that  at  all  stress  levels,  the  failure  times  are  governed  by  dis- 
tributions which  are  members  of  the  same  parametric  family,  such  as  an 
exponential  or  a Weibuil.  These  assumptions,  though  appealing  from  a 
statistical  point  of  view,  may  be  unreasonable  in  many  practical  situations. 
Of  particular  concern  is  the  last  assumption.  The  different  stress  levels 
may  have  different  effects  on  the  "failure  mechanisms",  and  thus  from  an 
engineering  point  of  view,  it  may  be  more  realistic  to  allow  for  different 
forms  of  the  failure  distribution  at  the  different  stress  levels. 

In  a recent  paper,  Proschan  and  Singpurwalla  (1978)  (hereafter 
referred  to  as  PS(l978))  have  proposed  a new  approach  for  making  inferences 
from  accelerated  life  tests  which  requires  neither  distributional  assump- 
tions nor  the  specification  of  a time  transformation  function.  Rather,  their 
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approach  is  Bayesian,  and  is  motivated  by  procedures  actually  used  in 
practice  by  reliability  engineers.  The  Bayesian  point  of  view  enables 
a user  to  incorporate  some  a priori  information  which  is  available  in 
accelerated  life  tests. 

In  Section  2 we  shall  review  this  pragmatic  Bayesian  approach  to 
accelerated  life  testing  by  PS (1978),  and  present  it.  in  a manner  which 
makes  it  easy  fc0  understand.  In  Section  3 we  shall  demonstrate  the 
usefulness  of  this  approach  by  applying  it  to  some  real  life  data  on 
accelerated  life  tests  presented  in  Nelson  (1970). 

Our  objective  in  writing  this  paper  is  to  make  this  new  approach 
to  accelerated  life  testing  accessible  to  the  reliability  engineer, 
and  to  demonstrate  its  usefulness  by  applying  it  to  a realistic  situation. 

2*  A-Pragmatic  Bayesian  Approach  to  Accelerated  Life  Testing. 

We  shall  .denote  the  k accelerated  stresses  (environments)  by 

EX,E2, . . . ,Efc,  and  the  normal  or  use  conditions  stress  (environment) 
by  E^.  Let 

E1  >E2>  \ >\  , (2.1) 

where,  E±  > E denotes  the  fact  that  E,  is  more  severe  than  E . 

Let  Fj  be  the  failure  distribution  of  the  items  tested  under 
environment  E..,  and  let  (x)  denote  the  failure  rate  of  F.  at 

ti 

time  point  x > 0. 

Because  of  condition  (2.1)  it  is  logical  from  a physical  point  of 
view  to  assume  that  for  any  x > 0, 
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(2.2) 


\ (x)  > Xg(x)  > >\k(x)  >\u(x)  . 

Using  failure  data  obtained  under  E-^E^ . . . ,Ek,  we  would  like  to  obtain 
Vj(x)>  an  estimate  of  X..  (x),  j = 1,2,..., k,  such  that  for  some 
0 < L < oo,  and  all  x e [0,L],- 

^ St  /V  St  St  A 

\(x)  > kg(x)  > >\k(x)  . (2.3) 

st 

The  notation  X > Y denotes  the  fact  that  X is  stochastically 
larger  than  Yj  that  is 

P[X  > x]  > P[Y  > x]  for  all  x. 

In  order  to  obtain  estimators  X.  (x),  j = 1,2,..., k,  which 
satisfy  (2.3),  we  shall  use  a Bayesian  approach.  Under  this  approach, 
condition  (2.2)  is  incorporated  as  a prior  assumption.  However,  we 
shall  first  define  the  "average  failure  rate"  and  discuss  its  Bayesian 
estimation  unconstrained  by  (2.2). 

2.1.  A Bayesian  Estimation  of  the  Average  Failure  Rate. 

Let  3^j(t)  be  the  number  of  items  undergoing  life  test  at  time  t, 
under  environment  E . Divide  the  time  interval  [0,L]  into  sub-intervals 
of  length  h > 0,  where  h is  chosen  to  make  L a multiple  of  h.  Thus, 
the  total  number  of  sub-intervals  is  L/h.. 

For  convenience,  we  shall  introduce  the  following  notation: 
ti  = the  time  point  (i-l)h,  i = 1,2,. . . , (L/h) 

[t^ti+h)  = the  i^*1  time  interval 
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N.  . — the  number  of  items  exposed  to  the  environment 


E.  at  time  t. 

0 l 

X.  . - the  failure  rate  of  F.  at  t. 

Oj3-  Oi 

± = ttie  number  of  failures  in  the  ith  time  interval 

under  environment  E . 

3 

P0,i  = the  Probabili^y  of  failure  of  a unit  in  the  ith 
time  interval  under  environment  E.. 


Let 


p.  . 

JL£- 


i-1 

- £ 

i=i 


-3}Z 


then  as  the  conditional  probability  that  an  item  which  has  survived 

to  time  t^  will  fail  by  time  tj[+h.  We  can  also  interpret  p*  . as 

j*1 

the  average  failure  rate  over  the  interval  Ft. .t.+h). 

i7  i ' 

If  we  assume  that  there  are  no  withdrawals,  removals,  or  censoring, 

then 


* = N.  . 

l+l  o,i 


i — 1,2, . » . , (l/h ) . 


Our  Bayesian  analysis  involves  assigning  prior  distributions  to 
the  average  failure  rates.  Suppose  that  the  p*  . , i = 1,2, . . . , (L/'h), 
are  a priori  independent  beta  random  variables  with  prior  parameters 
a and  P so  that  their  marginal  densities  are 

ti 


t(v*  ) _ r(^o}  / * f-1  , * xpo_1 

f(p0,i}  r(a;r (p.)  (po,i}  ^i,^ 


k 


•®ienJ'  given  the  N.  .’s  and  the  x.  . 's,  the  posterior  density  of 

O 9 A J * 1 v 

(pL^PJ,2-"‘^J(L/h))  is  [cf*  I’S(1978)1. 


(L/h) 

7T 


r(a+e.+N  ) 


1 ' r(a+x  ,)r(p.+N.  <~x,  J 


1=1 


o,i 


0 0,1  11' 


CX+x  —1 

(Pji)  ^ (l-P*.) 

0,1 


P.+N.  -■ 

0 0,1 


•x . . -1 

0,1 


(2.4) 


2*2.  A Bayesian  ijiStimation  of  the  Ordered  Average  Failure  Rates. 

In  the  Bayesian  context,  condition  (2.2)  leads  us  to  the  requirement 
that  for  every  fixed  value  of  i. 


* 


-l,i 


for  j = 2,3, . . . ,k  . 


Thus,  the  prior  distributions  of  p.  » . and  p*  . will  have  to  be 

0 **i,i  0,1 

chosen  so  as  to  ensure  that 


P^pJ-l,i  ~ > p£pJ,i  > p]  for  all  P,  0 < p < 1 . (2.5) 


One  way  of  achieving  the  above  condition  is  to  require  [see  FS(1978), 
Appendix  A]  that 

P.  > P.  - for  j = 2,3,..., k . 

To  be  assured  that  condition  (2.5)  is  also  satisfied  with  respect 
to  the  posterior  distributions  of  p*  and  p?  . , it  is  sufficient 

Qji  — 

that  for  every  fixed  value  of  i, 
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and 


x . . < x . , . 
0 , 3 —1 ) i 


(2.6) 


P.+N.  . -X.  . > , +TJ  -X 

3 J,i  3,1  — j-l  3-1, i 


for  3 = 2,3,.. .,k. 

.The  first  part  of  condition  (2.6)  states  that  the  number  of 

failures  in  the  interval  [t^^-fh)  under  environment  E„  must  not 

J 

be  greater  than  the  number  of  failures  in  the  same  interval  under 
environemnt  for  all  values  of  3*.  Further-more,  as  is  discussed 

in  PS  (1978)  a reasonable  strategy  is  to  have  more  items  initially  on 
test  under  the  more  severe  environment  than  under  the  environment 

Ej,  so  that  ^ for  all  3. 

Since  N^i  = N3,i-i~X3^i._i>  ^e  second  part  of  condition  (2.6) 
can  be  written  as 


VfW+i  s y»3,1+1 


for  3 = 2,3,..., k 


Thus  for  every  fixed  value  of  i,  the  number  surviving  at  the  start  of 
s *fc 

the  (i+l)  interval  plus  the  prior  parameter  for  the  environment 

E3-l  must  be  not  greater  than  the  corresponding  sum  for  the  environment 
E^..  Since  the  number  of  failures  in  a particular  interval  is  a function 
of  the  severity  of  the  environment  and  the  number  on  test,  and  since 


a reasonable  strategy  is  to  choose 


6 


Thus,,  the  prior  parameters  S nr,^  r , . 

pj„i  EB.  Pj  are  indicative  of  the 

£g^-^e^J^rfcy_of  the  environmental  conditions  e and  E . . 

Because  of  condition  (2.1): 


*1<*2<  <K  , 


vdth  the  values  of  tte  fys  being  indicative  of  the  severity  of  the  B 

3 

If  the  prior  parameters  Pj,  the  failure  data  Xj  and  the  N. 

are  such  that  condition  (2.6)  is  satisfied  for  eveiy  fired  value  of  i^ 
then  the  stochastic  ordering  condition  (2.5)  will  be  automatically 
satisfied  »ith  respect  to  the  posterior  distribution  of  n*  w* 

f this  is  not  the  case,  then  we  shall  "pool”  the  adjacent  violators 

using  the  pooling  procedure  described  below;  the  purpose  is  to  eliminate 
violations  of  (2.6). 


2,5“  Pie  Pooling  of  Adjacent  Violators. 

Hie  procedure  for  pooling  adjacent  violators  described  here  is 
commonly  used  in  isotonic  regression;  see  Barlow,  Bartholomew,  Bremner 
and  Brunk  (1972). 

Consider  the  time  interval  [(i-l)h,ih);  by  condition  (2.6)  we 
require 


Xl,i  - -•  > 


-x2,i-"'  > xi_i  ,*  > x . . > • . . > v 

* J-l^i  ~ j,x  - - Tc,i 


and 
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Pi+Vxz,i  ?-  P2+Vx2,i  < 


< P.  , +N.  , .-X.  , . 

- 3-1  j-l,i  j-l,i 


2 W-xj,i  s •••  < w-*ki 


If  a reversal  occurs,  that  is,  if  either 


X.  - . < X.  . 
3«jLjX  0*3- 


or  if 


Pi-1+NT  1 4“**  1 A > P,+W.  .-X.  . , 

3 -*•  3 r*-j3-  J "*_L « X j * 

then  we  pool  the  violators  and  replace  them  as  shown  below. 

Replace  both  x..^.  and  ^ by  |{x._^i+x.^ . ) mS 

b(&i  i +£ . +N  , +N.  . -X,  ,-x.  . ). 

d 0-1  3 J~l,i  J>i  j-l,i  j,i' 

We  now  test  the  new  sequences  to  see  if  they  are  properly  ordered. 

If  not,  we  replace  the  adjacent  violators  by  their  appropriate  averages 
Thus,  if 

2^xi-l  i+x-i  *1  ) = k(xn  i ,+x.  .)  < x.^,  . , 

^ J x,i  j,i  d 0-l^x  o+l,i  * 

then  we  replace  each  of  the  three  by  the  average 

J(xj-l,i+xj,i+xj+l,i)  * 


The  same  procedure  is  used  for  the  (p.+N.  .-x.  . )'S. 

0 <3>x  J,i 
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We  follow  the  above  procedure  for  all  time  intervals  and  continue 

until  all  reversals  are  eliminated.  Note  that  excessive  pooling  will 

occur  if  the  relationship  (2.2)  does  not  actually  hold,  or  if  the 

environmental  conditions  are  too  similar  to  each  other. 

Assuming  a squared  error  loss,  before  pooling,  the  Bayes  estimator 

of  p*  is  [see  PS(1978)] 

3 A 


'St* 


a+x.  . 

a+e.+N.  . 

3 3#i 


(2.7) 


If  we  have  done  some  pooling,  then  the  x.  .'s  and  the  ({J.+IT.  . )'s 

3^1  3 3,1 

are  replaced  by  their  appropriate  pooled  averages. 

2.4,  A Made!  for  Extrapolation  to.  Use  Condition  Stress. 

We  shall  use  the  to  estimate  P*  rP*>2,...,P*;  (L/h), 

the  average  failure  rates  under  the  use  conditions  environment  E . 

u 

In  the  absence  of  any  physical  or  engineering  knowledge  about  the 
relationship  between  the  average  failure  rate  and  the  corresponding  stress, 
we  shall  postulate  the  following  simple  but  reasonable  relationship  among 
the  average  failure  rates. 

For  each  value  of  i,  i = 1,2. . , . , (h/h). 


* 


PC 


-k,i  0 l^k-1, i 


. + 


* 


+ Wk-lPl,i  ' 


(2.8) 


where  wq,w1^  * * * ^k-l  are  un-^nown  constants. 

The  above  relationship  states  that  the  average  failure  rate  over 
a particular  time  interval  under  the  environment  E is  a weighted  sum 

xC 

of  the  average  failure  rates  over  the  same  time  interval  under  the 


environments  E^  ^,E^  2, . . . ,E^ 
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_ , ^ A A 

LCo  * * * ,wk-l  ’ ”^e  Teast-squares  estimators  of 

wq>w1>  * * * >wk*  these  can  be  obtained  routinely  from  equation  (2.8), 
but  with  p*  in  place  of  p*  throughout.  Thus  for  i = 1,2, . . . , (i/h), 
we  have 


P*  . 
u,i 


A A Ji. 


+ 


k-l*2,i 


as  the  estimators  of  the  average  failure  rate  under  E . As  a conse 
quence,  we  also  have 


Pu, (L/h)+l  “ 1 


(L/h) 

I 

i=l 


An  estimator  of  F^(t),  the  probability  of  an  item  surviving  to 
time  t under  E^,  the  use  conditions  stress,  is 


3.  An  Illustrative  Real  Life  Example. 

In  this  section  we  shall  apply  the  methodology  discussed  in  the 
previous  section  to  some  accelerated  life  test  data  given  by  Nelson 
(1970).  These  data  represent  the  times  to  breakdown  of  an.  insulating 
fluid  subjected  to  elevated  voltage  stress  levels.  For  convenience,  we 
shall  consider  here  only  4 accelerated  voltage  levels  - 36,  34,  32,  and 
30  kilovolts  (KV's).  The  failure  times  (in  minutes)  under  the  various 
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stress  levels  are  given  in  Table  1.  Kelson* s original  data  correspond 
to  7 different  stress  levels,  but  some  of  these  contain  very  few  failure 
times  and  are  therefore  omitted. 

We  shall  assume  that  the  use  conditions  stress  is  28  KV,  anri  apply 
our  approach  to  estimate  the  failure  distribution  of  the  breakdown  times 
at  this  stress. 

Following  the  notation  of  Section  2 , we  shall  choose  L to  be  100 

minutes , and  h to  be  1 minute,  making  the  total  number  of  time 

intervals  equal  to  100.  For  our  prior  parameters  we  shall,  following 

the  discussion  in  Sections  2.1  and  2.2,  choose  a = 1,  with  p »i,  p = 2, 

and  P^  = 12.  The  above  choice  of  our  prior  parameters  p } p f 

X 2 

Py  and  P^  is  motivated  by  an  inspection  of  the  failure  times  in 

Table  1.  Observe  that  even  though  the  stress  levels  decrease  successively 

by  only  2 KV*s,  the  corresponding  change  in  the  failure  times  from  one 

stress  level  to  the  next  lower  appears  to  be  more  drastic.  This  may  be 

due  to  a change  in  the  failure  mechanism  as  we  go  to  the  lower  stresses. 

Note  that  ^ = 15,  K2  = 19,  ^ =15,  and  = 11. 

After  computing  the  x . *s  and  the  (P.+N,  .-x.  „)'s  for 

3 3,1 

3 — 1,2,3,^,  i = 1,2,..., 100,  and  after  pooling  the  adjacent  violators 

as  discussed  in  Section  2.3,  we  obtain  the  Bayes  estimators  p*  . using 

3,2. 

Equation  (2.7).  These  values  are  given  in  Table  2,  wherein,  for  convenience, 

we  have  limited  ourselves  to  time  intervals  of  width  5.  Observe  that 

the  pooling  scheme  ensures  that  for  each  time  interval  the  n*  * s 

3,i 

decrease  with  increasing  stress  level.  For  example,  p*  = .3333 

20  ^ 

^1,20  = *125°>  £3,20  = ,071^  ^ i*,20  - *0454' 
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Our  next  step  is  to  use  the  entries  in  Table  2 to  obtain  estimates 
of  the  w’ s given  in  Equation  (2.8).  If  we  designate  the  30  KV  level 
as  the  level  k(=4-)  of  Equation  (2.8),  and  the  32,  3^,  and  the  36  KV 
levels  as  k-1,  k-2,  and  k-3,  respectively,  then,  the  least  squares 
estimation  of  the  w' s gives  us  for  i = 1,2, ...,100, 

± - .2360  + .I886ij*  + . 33530P* ±-+  .09857?^  * 

'Si. 

¥e  shall  now  use  the  following  equation  to  estimate  p •., 

"U.  1 

i = 1,2,..., 100,  where  u(=5)  denotes  the  use  condition  stress 
level  of  28  KV: 

P5  ± = .2360  + .l886p^i  + .13538p^i  + .09857?*^  * 

An  estimate  of  the  survival  probability  at  use  conditions  voltage 
of  28  KV  is  . given  by 


Hie  values  of  F^(t),  t = ^,10,  • * • >100,  are  given  in  Table  3> 
and  plotted  in  Figure  1. 
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Table  1. 


Times  to  Breakdown  of  an  Insulating  Fluid 
(in  Minutes ) Under  Various  Values  of  the 
Stress  (in  Kilovolts). 


36kv 

34KV 

.35 

•19 

• 59 

.78 

.96 

.96 

• 99 

1.31 

1.69 

2.78 

1.97 

3.16 

2.07 

4.15 

2.58 

4.67 

2.71 

4.85 

2.90 

6.50 

3.67 

7.35 

3-99 

8.01 

5*35 

8.27 

13.77 

12.06 

25.50 

31.75 

32.52 

33.91 

36.71 

72.89 

32KV 

3QKV 

.27 

7.74 

.40 

17.05 

.69 

20.46 

.79 

21.02 

2.75 

22.66 

3.91 

43.40 

9.88 

47.30 

13.95 

139.07 

15.93 

144.12 

27.80 

175.88 

53.24 

194.90 

82.85 

89.29 

100.58 

215.10 
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Table  2 


• .Estimated  Values  p_.  . of  the  Average 

...  D s i . 

Failure  Rates  Under  Various  Values  of 
the -Stress  (in  Kilovolts). 


Stress 


Time  Interval 

36kv 

34kv 

32KV 

30KV 

. 5 

.3846 

.1724 

.0588 

.o4i7 

10 

.3077 

.1429 

.0816 

.0435 

15 

.3333 

,1250 

.0667 

.0435 

20 

.3333 

.1250 

.0714 

.0454 

25 

.3333 

.1250 

.0714 

.0526 

30 

. 5000 

.1250 

.0769 

.0526 

35 

.5000 

.2000 

.0769 

.0526 

4o 

.5000 

.2500 

.0769 

.0526 

45 

.5000 

.2500 

.0769 

.0556 

50 

.5000 

.2500 

.0769 

.0588 

55 

.5000 

.2500 

.0833 

.0588 

6o 

.5000 

.2500 

.0833 

.0588 

65 

.5000 

.2500 

.0833 

.0588 

7° 

. 5000 

.2500 

.0833 

.0588 

75 

.5000 

.3333 

.0833 

.0588 

8o 

. 5000 

.3333 

.0833 

.0588 

85 

.5000 

.3333 

.0909 

.0588 

90 

.5714 

.4000 

.1290 

.0588 

95 

.5000 

.3333 

.1000 

.0588 

100 

.5000 

.3333 

.1000 

.0588 
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A 

Table  3.  Values  of  F(t)  at  the  28  KV  level,. 


A 


t 

F(t) 

5 

.803 

10 

.650 

15 

.525 

20 

.424 

25 

.343 

30 

.277 

35 

.222 

40 

.177 

45 

.141 

50 

.112 

55 

.090 

6o  ‘ 

.071 

65 

.057 

70 

.045 

75 

.036 

80 

.028 

85 

.022 

90 

.017 

95 

t 

0 

H 

■{=- 

100 

.011 

*5 


Survival  Probability 


Time, 
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